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ABSTRACT 

The deposition of mechanical feedback from a supermassive black hole (SMBH) in an active galactic 
nucleus (AGN) into the surrounding galaxy occurs via broad-line winds which must carry mass and 
radial momentum as well as energy. The effect can be summarized by the dimensionless parameter 

77 = Moutf/-^acc = ZewC^/i'w where Ew (= Evj / (MaccC^)) is the efficiency by which accreted matter 
is turned into wind energy in the disc surrounding the central SMBH. The outflowing mass and 
momentum are proportional to rj, and many prior treatments have essentially assumed that 77 = 0. 
We perform one- and two-dimensional simulations and find that the growth of the central SMBH is 
very sensitive to the inclusion of the mass and momentum driving but is insensitive to the assumed 
mechanical efficiency. For example in representative calculations, the omission of momentum and mass 
feedback leads to an hundred fold increase in the mass of the SMBH to over 10^^ Mq. When allowance 
is made for momentum driving, the final SMBH mass is much lower and the wind efhciencies which 
lead to the most observationally acceptable results are relatively low with ew S 10^^. 
Subject headings: accretion, accretion discs - black hole physics - galaxies: active - galaxies: nuclei - 
galaxies: starburst - quasars: general 



1. INTRODUCTION 

Feedback from active galactic nuclei (AGNs) at the 
centers of galaxies is believed to have a significant effect 
on the evolution of those galaxies. However, the precise 
physical mechanisms by which this feedback occurs are 
greatly uncertain — perhaps more so than is commonly 
acknowledged. While much path-breaking and insight- 
ful work has been done, it is also true that some of the 
most basic requirements, such as the necessity that mass, 
energy, and momentum be conserved, have not been im- 
posed in several of the popular treatments of this subject. 
And the inclusion of the presently known and observed 
feedback processes is often treated selectively. The pur- 
pose of this paper is to attempt to lay out the physical 
framework for discussing the issues and to provide illus- 
trative examples of the results obtained primarily from 
one-dimensional computations that include or exclude 
specific processes. We also include a treatment of the 
two-dimensional, axisymmetric case, presented in less de- 
tail, to show how the qualitative features carry over to 
this more realistic case. Definitive solutions are beyond 
present art in this field, so the focus will be on the qual- 
itative features of the physical solutions rather than the 
detailed comparison with observations. 

In outline, there are three phases to the overall phe- 
nomenon: 1) the provision of fuel for the central super- 
massive black hole (hereafter SMBH); 2) the generation 
of the outflowing stream of energy, mass and momentum 
from the vicinity of the SMBH; and 3) the absorption and 
transmission of this energy, mass and momentum by the 
ambient gas in the galaxy and the subsequent reactions 
of the ambient gas to that input. 

I. The fueling is generally believed to be via in- 
falling gas, and typically, two origins for that gas have 
been proposed; at high redshift ambient gas in discs 



liberated by the non-axisymmetric forces released dur- 
ing m ergers is certainly important ([Barnes fc Hernquistl 
Il991t ) , while at lower redshift merg ers fail and proba- 
bly do not fuel AGNs ()Li et al.ll200l( ). but the processed 
gas released via normal stellar evolution p rovides an am 
pie source ( Mathews 1983; ShuU 1983; Ci otti et al.l ll99I 
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iPadovanlfc Matteuccii iI993c iCiotti fc Ostrikeii 
hereafter CO07). A primary clue as to which of these 
sources dominates in a specific case is provided by de- 
tails of the metallicity distribution, since the reprocessed 
gas probably has super-solar metal abundance and will 
also show signs of stellar evolution such as higher ni- 
trogen or S-process abundances. The clue to fueling by 
infalling globular clusters might be relatively low abun- 
dances of elements made in SN I such as Fe. In almost 
all treatments a central disc mediates between the in- 
flowing material and the SMBH. Other sources, such as 
small stellar systems dragged in by dynamical friction 
have been considered from time to time. These stars, 
or oth ers added to the central regions via loss-cone pro- 
cesses fBe gelman. Blandford fc"Reesl lT98^ can be shred- 
ded during tidal interactions with the central SMBH or 
by collisions with one another or with a central disc — 
the debris collecting in the disc and fe eding the cen- 
tral SMBH via conventional m echanisms (|Ostrikeiill983l : 
IDai. Fuerst fc Blandfordl[2010h . 

2. The outflows fall into three categories. The sig- 
nature of AGNs is, of course, the enormous electro- 
magnetic, luminous output, with major contributions 
from the IR bands to the gamma ray region. The 
bulk of the flux is typically in the "UV bump" and the 
flux from this region thus doi ninates for the "momen- 
tum driven winds" (e.g., see [P roga. Stone fc K allman 
20M lKh^[2003L iProga fc KaTh^an,2004: DeBuhr et al 



20091) . The region where this driving occurs is fairly 
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close to the quasar (50i?BH ^ f < 500i?BH, where 
i?BH is the SMBH Schwarzschild radius). However, 
the moderately hard X-rays determine the average pho- 
ton energy: (hv) = h J L^vdv/ ^ L^Av when integrated 
over the spectrum. This region of the spectrum domi- 
nates the photon heating, causing the heated gas to ap- 
proach the mean photon energy: l.SfcTx = {hv) with 
Tx - 2 X lO^K (ISazonov. Ostriker fc Sun^^ [Mil . 
The resultant heating occurs over an extended range 
of radii: 1 00 pc< r <3 kpc (jCiotti. Ostriker fc Progal 
120091 120101 hereafter Papers I and III ), and it can be 
significant for r > 0.1 pc (jProgal 12007). It can effi- 
ciently drive outflows as shown in a series of papers by 
Ciotti, Ostriker and collaborators (cf.COO? and refer- 
ences therein). For electromagnetic output there is, of 
course, no rest-mass component. The total energy emit- 
ted in this form ha s been established fairly accurately via 
the iSoltanI (|1982[ ) argu ment to be = ep ,]AM„,r.r.c^ 
with erad 0.1 - 0.15 (|Yu fc Tremainel[200a) . The mo- 
mentum output, of course, is Ap = AE/c. In opti- 
cally thick cases (r ^ 1), the total mome ntum absorbed 
by th e f luid can appro ach Ap = tAE/c (jPeBuhr et al.l 
|2009[) . iSilk fc Nusse^ (, 201 0) also consider the impor- 
tance of radiative momentum driven winds on galactic 
and cluster scales but limit the input to L/c, which can 
be con siderably less than allow ed in the optically thic k 
case bv lCiotti fc Ostrikeri ()2007[ ) orlDeBuhr^etalJ (|2009[ ). 

Next, let us turn to mechanical output. Both broad- 
and narrow-line regions inject mass, energy and momen- 
tum into the surrounding gas, with the broad-line winds 
probably dominant. Since these are material flows with 
velocity in the vicinity of the SMBH, v^, the mass out- 
flow can be considerable. If we let the inflowing and 
outflowing mass rates be (Afi„f, Afoutf), then conserva- 
tion of mass, energy and momentum can be summed up 
with the following simple equations: 



Mi. 



OUtf ; 



(1) 



where Mace is the mass rate actually accreted by the 
SMBH, and 



-E'w = ^MoutfWw 
Pw = Moutffw, 



(2a) 
(2b) 
(3) 



are the wind energy and momentum, respectively. We 
have oversimplifled matters by allowing only one wind 
velocity, when in fact Equation [^a] requires (w^) and 
Equation [3] requires (ww)- Also, it is important to specify 
exactly where and when the quantities in equations [IH3] 
are to be measured. In the conventional treatment of this 
subject, the SMBH is surrounded by a disc or torus to 
which matter has fallen from larger radius. Then, plac- 
ing a sphere around this disc or torus (at r ~ 1 pc), 
the instantaneous spherically averaged infall through the 
sphere is Mini{t) and the spherically averaged outflow 

is Moutf(^)- The difference will be accreted onto the 
SMBH unless driven out in disc originating winds; the 
latter of course contributes to Afoutf and so the remain- 
der Mint — Moutf will be accreted. Two further compli- 
cations are allowed for in some detailed treatments: (i) 



the actual instantaneous value of Afacc is a time-lagged 
convolution of the quantity in Equation [T] since a finite 
time elapses as material is transported through the disc 
to the central SMBH and (ii) star formation may (in fact 
frequently will) occur in the disc, removing mass that 
would otherwise have accreted onto the SMBH. Both of 
these complications are allowed for in CO07 and other 
work, neither is of dominant importance. 

N ow, defining the dimensionless ratio from Eauation l2al 
and [2bl to be. 



M 



outf 



2e„ 



we can now rewrite Equations [TH3] as 



Macc = Mi„f 

Moutf = A/inf 
1 



1 



1+77' 

l + ry' 
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CwA/infC 



l-t-ry 



1+77 



(4) 

(5a) 
(5b) 
(5c) 
(5d) 



These Equations, [S^-d, are, in fact, the ones that most 
authors have adopted who treat AGN feedback as a uni- 
fied process comprising both infall and outfiow. How- 
ever, they typically adopt = 0, implicitly assuming 
Ww — >■ 00, so that Moutf and Pw are neglected and the 
two terms that are included, E^ and Mace, may be 
overestimated. If it eventuated that ry really is a very 
small number, then not much error would be induced 
and one would be justified in neglecting the out-flowing 
mass and momentum and in setting E^ ~ ewMjnfC^, as 
most authors assume. If we adopt for the efficiency of 
generating mechanical ener gy the value Cw = 5 x 10~p, 
as done bv iSpringel. Di Matteo. fc Hernguistl (j2005D : 
[Johansson. Naab fc Burkert I (I2009D. hereaf t er SDM H05 
and JNB09 respectively, iMcCaxthiLeLaLl dMS) and 



other authors, and we take 



10'' km/s (z;w,io) 



(|Moe et al.l 120091) , then we have from Equation HI r\ 
9w~^jg and all of the neglected effects may in fact be 
dominant; the bulk of the inflowing mass may be ejected 
in a broad-line disc wind, and the mass and momen- 
tum input deposited in the ambient gas may dominate 
over the energy input, which may be largely radiated 
away. Papers I and Paper III do include these effects, 
but do not spell out their signiflcance. The principal 
purpose of the present paper is to do just that — to show, 
in specially simple one- and two-dimensional calcula- 
tions the effects of including or excluding mass, energy 
and momentum conservation when 77 > 0. In addition 
to the papers referred to above which attempt to com- 
pute both the infall to the central SMBH and the out- 
flow from it in a unifled fashion, there are many others 
which postulate a central source and then, after esti- 
mating the mass, momentum and energy flowing at of 
that source (and some angular and temporal distribu- 
tion thereof) do effectively compute the effects of that 
injection of energy, mass and momentum onto the sur- 
rounding fluid. Space does not permit a comprehensive 
description of this related subject of research, but impor- 
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tant p a per s includ e the fol l owing : cf. iMetzler fc EvrardI 
([T99I: ISte rnbcrg fc Sokeri (12003): iFabian et al.' ' 
Reeves et al.i (i2009D: lArieli. Rephaeh fc NormanI ()2010[ ): 



Gaspari et all (|2010( ) 

The wind efficiency, ew, is not known very well- 
neither from observations nor from detailed physi 
cal simulations. But the best estimates from ei 
ther of these sources might be in th e range 1 



IQ- 



> e 



> 3 X 10"^ (iProga. Stone fc Kallma n I2QOOI: 

Proga fc Kallman"2004'; 'Krongold et al."2007i; iStoll et alJ 
2009: Kurosawa, Proga & Nagaminc 2009!) , a factor of 5 
to 17 smaller than the commonly adopted values and in 



a range where 77 < 1 if t > w.io ?a 1. A specific example 
may be useful. iMoe et al.l (|2009D study the quasar SDSS 
J0838+2955. They find a mechanical energy output of 
4.5 X 10''^ erg/s, a mass outflow rate 10 times the accre- 
tion rate and a mechanical efficiency of 1 x 10~^, and they 
quote other observational studies which indicate similar 
numbers. From analyses of the ionization parameters in 
the broad-line winds, estimates of the radial extent of the 
winds can be made; the above paper, and those quoted 
within indicate radii measured in kpc — consistent with 
the one-dimensional numerical work in Paper III. 

As shown in Papers I and III, an additional important 
question asks "what fraction of the sky is covered with 
the broad-line winds?" Again two approaches are possi- 
ble. Empirically, on the order of 20-25% of bright quasars 
show broad-line winds; this translates to ~ tt steradians 
or 7r/2 steradians in each conical outflow, if we assume 
that the wind is emitted symmetrically above and below 
the inner AGN disc. On the theory side , the r adiation 
driven winds found by IProga fc KallmanI (pOOl . via de- 
tailed hydro radiation-transfer calculations, cover ~ tt 
steradians, roughly consistent with the observational es- 
timates. 

Finally, let us turn to the narrow jets, the outflow 
observed from AGN in "radio mode", when the elec- 
tromagnetic luminosity is considerably below the Ed- 
dington limit. M87 is an excellent nearby example of 
such a system. These are standard FRI radio sources. 
Here the jets are quite narrow and appear to be com- 
prised primarily of a relativistic fluid. The same type 
of calculation as presented in the last section would 
indicate that the out-flowing mass is of negligible im- 
portance and the energy output greatly dominates over 
the momentum output. The total energy output from 
these phases is considerable, but the accretion rates are 
thought to be low in these phases so the efficiencies of 
ener gy generation may be very h igh (cf. for a computa- 
tion jMcffinnexEGimmiillQOl) . Since so little mass is 
accreted in radio mode, the Soltan argument cannot be 
used to empirically estimate efficiencies, but, from the 
observational estimates of the energies available in the 
giant radio lobes, it may be that the AGN emits in ra- 
dio mode considerably more energy than it does in wind 
mode. However the deposition from the intense but ex- 
tremely narrow streams appears to be inefficient, and 
the jet drills through the gas in the surrounding galaxy, 
dumping most the energy into the intergalactic medium. 
Thus, while it may act as the dominant feedback mecha- 
nism for the IGM (and we will return to this in a subse- 
quent paper) , it is of lesser importance than the radiative 
or wind components in heating and driving out the am- 



bient gas from within a galaxy. 

3. The interactions between the out-flowing energy, 
mass and momentum with the ambient fluids are com- 
plex and are just beginning to be studied with the needed 
detail. We focus here on the relatively gas poor ellip- 
tical systems, since it is in these that the bulk of the 
mass in SMBHs is found. The radiative interactions are 
perhaps easiest to describe. Since the mechanical mo- 
mentum is conserved and cannot be radiated away, it 
can be a dominant effect. The minimum level of inter- 
action is provided by electron scattering and, since the 
most luminous quasars are found to be clustered near the 
Eddington luminosity limit (at which level the momen- 
tum absorbed by electron scattering balances the gravi- 
tational force on the fluid from the central SMBH), we 
know that this effect is significant in many cases. Ab- 
sorption of the out- flowing radiation will not, in gen- 
eral, reduce this effect, since typically the radiation is 
simply re-emitted in another band and electron scatter- 
ing opacity is wavelength independent until the Klein- 
Nishina limit is reached at very high energies. In fact, in 
the optically thick limit, the radiation is transformed by 
dust absorption into the infra-red, but the effects in this 
case are even greater than in the simple case, since the 
scattering opacity of the dust to infra-red is, per atom, 
larger (by roughly a factor of 5) than the electron scat- 
tering cross section. For the bright ULIRGs, which may 
contain both an active AGN and a brighter starburst, 
there will be a near balance between the inward grav- 
itational forces and th e outward radiative momentum 
trans fer on the dust (cf. [Thompson. Quataert fc Murray! 
120051 COOT). Under these circumstances, the inner sev- 
eral hundred parsecs of the galaxy are analogous in their 
equilibrium structure to a very massive star in so far as 
there is a nearly equilibrium balance between radiative 
and gravitational forces. 

The effects of heating from the AGN are, for quite dif- 
ferent reasons, also likely to be independent of absorp- 
tion (so long as it is not excessive, i.e. not Compton 
thick). Sazono vet al.! ([2005) present a simple analytical 
exploration of the effects and Paper I presents a more 
detailed one-dimensional treatment. The photons which 
dominate the heating process are in the moderately hard 
region (^50keV), and we know from X-ray absorption 
studies that AGN are typically optically thin to such ra- 
diation. Metal line resonance absorptions dominate the 
absorption unless the spectrum is extremely hard, and 
in those cases Compton absorption would be dominant. 
If we consider the issue on a per atom basis, all that 
matters is the heating per atom, which scales as (as- 
suming that the fluid is optically thin to hard X-rays), 
and the cooling rate per atom which scales as the density. 
Since the latter can also scale as or even falls off at a 
steeper rate in some circumstances, the heating can bal- 
ance or exceed cooling over an extended range of radii. If 
that happens, the gas temperature will rise towards the 
radiation temperature, Tx ~ 2 x lO"^ K. Then, since this 
exceeds the virial temperature in almost all galaxies, the 
heated gas, having thermal energy higher than its gravi- 
tational energy, will be accelerated outwards and tend to 
drive a wind into the surrounding fluid. Of course, since 
this will shut off the accretion flow and the fuel to the 
central source, the result will be a burst of energy output 
followed by much slower cooling of the shocked gas and 
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a repeated burst at a much later time. Thus, episodic 
accretion is expected. 

The mechanical energy input is more localized to the 
vicinity of the SMBH and would be efficient in "protect- 
ing" the SMBH from very high rates of accretion, except 
for one important caveat. It necessarily happens that 
such episodes of high rates of energy deposition will oc- 
cur when the central gas densities are high, and under 
such circumstances the gas will tend to radiate away the 
input energy unless f orbidden t o do so as has occurred 
in some calculations (jBooth fc Schayc 2009). This, as 
we shall see, makes the energy input rather inefficient in 
driving outflows and in protecting the SMBH from ex- 
cessive accreation. But the momentum input cannot be 
radiated away, and, as we shall see in the remainder of 
the paper, it is very efficient in limiting the infall and 
accretion onto the central SMBH. Mechanical input, via 
either thermal or momentum based mechanisms, will also 
tend to produces episodic accretion. 

The broad-line gas outflow must drive a strong 
shock into the ambient gas, and that, in turn, given 
standard physics, should accelerate charged parti- 
cles efficiently via a variant of th e first order Fermi 
process fcf. fB landforT fc Ostriked [TqtI iBell ITOTSl : 
iBlandford fcEichler. 19871 ). Then this relativistic fluid 
will further drive the outflow and, since thermal radia- 
tion is suppressed for this component, the conversion may 
somewhat enhance the effects of feedback. But, overall, 
this process simply transforms internal energy from one 
form to another and so, whereas it may be observation- 
ally quite significant, it will have a relatively small global 
effect. T wo recent p apers that h ave explored these pro- 
cess es arelFuiita et a l. (2007) and lJiang et al.l (|2010[) : see 
also ISironi fc So'crai;esl poioh . 

2. THE MODELS 

In this section, we summarize the main properties of 
the galactic models adopted in this study. A detailed 
description of the galaxy models and input physics is 
given in CO07, Papers I and HI. 

We study galaxy models characterized by the effec- 
tive radius of the galaxy stellar component i?c = 6.91 
kpc, an initial stellar mass Af» = 2.87 x 10^^ Mq, and 
central aperture velocity dispersion da = 260 km s~^. 
This represents approximate ly the typical galaxy 
which lYu fc Tremaind ()2002l ) find contain the bulk of 
the identified mass in SMBHs. The stellar mass distri- 
bution which i s emb edded in a dark matter halo is de- 
scribed by the iJaffd (fl983i ) model while the total mass 
density distribution follows a profile; all the rele- 
vant dy namical quantities need in the sim ulation are 
given in [Ciotti. Morganti fc de Zeeuwl (|2009 |). The ini- 
tial SMBH mass is Mr h = 2.87 x 10*^ M©, following the 
iMagorrian et al.l ()1998h relations (A/bh ^ lO~^Ah). The 
simulations are for an isolated, giant elliptical galaxy 
without the effect of the intracluster medium, as out- 
flow boundary conditions are set at the galaxy outskirts 
250 kpc), so that the interstellar medium (ISM) is 
provided by the recycled gas produced by stellar evolu- 
tion. The simulations starts at 2 Gyr, which corresponds 
to a redshift of z ^ 3.2 for the ACDM cosmology with 
Qm = 0.3, = I - ^rn ^ 0.7 and Hq = 70 km s'^ 
Mpc^^ and ends at 14 Gyr. 

The input physics of the model is fully described in 



Paper I. Here we recall the most important aspects. The 
instantaneous bolometric accretion luminosity is 



(6) 



and we adopt an advection-dominated accretion flow 
(ADAF)-like radiative efficiency as 



Am 
1 + Am' 



Af, 



(7) 



Edd 



where Msdd — ^Edd/eoc^j and A is a free parameter 
so that cem ~ e^A m for m <g: A~^ . W e flx A = 100 
in our simulations f Na ravan fc Yil [l994[) . and we adopt 
for the peak EM efficiency e g = 0-1 or 0. 2 consistent 
with estimates based on the iSoltanI (|1982[ ) argument. 
In the treatment of radiation feedback, we consider 
the radiation pressure as well as heating/cooling feed- 
back j including photoionization, Compton and line hcat- 
ing (iSazonov. O striker & Sunvacv 2004; Sazonov et al] 
l2005f ). In accordance with both observations and theo- 
retical expectation, the transformation of accreted mass 
to electromagnetic energy output declines dramatically 
at low accretion rates. 

In the mechanical feedback treatment, the fiducial in- 
stantaneous mechanical luminosity of the disk wind is 
calculated as 



La 



fwMaccC + eilC (1 - /re 



M, 



dh* 



(8) 



where Cw is the mechanical efficiency of the wind, and the 
second term represents the energy associated with the 
Type II supernova (SNII) explosions of the high-mass 
stars in the circumnuclear disk (see Paper I, Equation 
(20) for details). Here, Mdh* is the current mass in the 
disc in high mass [M > 9>Mq) and r,h is their typical life- 
time. In this work, we restrict attention to the commonly 
assumed case of a constant value of Ew (e.g. SDMH05), 
which corresponds to Type A models in Papers I and HI. 
Physically, a fixed mechanical efficiency implies that the 
mass accreted by the central SMBH has a fixed relation 
to the mechanical energy flowing out of the central re- 
gions. We here neglect the jet effects, which are expected 
to be effective only in the low-luminosity, hot accretion 
phases at late-time evolution. The reference models (AO 
and Al) from Paper III study the evolution of gas and 
the mechanical feedback from SMBHs, and solve Eule- 
rian equations of hydrodynamics with mass, energy, and 
momentum sources (see Paper I). In order to study the 
effect of each physical process, i.e., mass, energy, and 
momentum feedback, we build several models which ne- 
glect one or two of physical terms. We discuss the details 
of each model and comparison of them in the following 
section. 

3. EXPLORING ONE-DIMENSIONAL MODELS 

The model properties and results are given in Ta- 
ble 1. The mechanical efficiencies Ew are given in Col- 
umn 3 and the corresponding values of = 2t^c^ /v'^ 
are given in Column 4. We devote Column 5-9 to 
present {z = 0) model properties. First, for some 
models indicated with the symbol ^ in Column 5, we 
distribute the mechanical feedback energy (momentum 
and mass) only at the lower boundary of the grids to 
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TABLE 1 

Summary of One-dimensional Model Properties 











Bottom 


Feedback 












Model 






Layer'' 


Radiation 


Energy 


Momentum 


Mass 


log AA/bh 


log Lx 


log Z|» = 


log Aiyw 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 




1 


AO ^ 


5 X lO^-' 


9 


X 




V 


V 


v/ 


6.72 


37.11 


-7.98 


58.67 


2 


MAO " 


5 X lO^-' 


9 


X 


X 


V 


V 


v/ 


6.72 


37.11 


-7.98 


58.67 


3 


EPMO-R 


5 X 10^^ 


9 


V 


V 


V 


V 


V 


7.07 


39.75 


-10.58 


59.02 


4 


EPMO 


5 X 10^-'' 


9 


V 


X 


V 


V 


V 


7.13 


37.84 


-10.58 


59.08 


5 


PMO 


5 X lO^-' 


9 


V 


X 


X 


V 


V 


7.13 


37.84 


-10.58 


59.08 


Q 


EO 


5 X 10^^ 


9 


/ 

V 


X 


V 


X 


X 


10.32 


41.31 


-3.61 


62. 2T 


7 


Al = 


2.5 X 10"* 


0.45 


X 


V 


V 


V 


V 


7.38 


36.36 


-6.72 


58.02 


8 


MAI ' 


2.5 X 10"^ 


0.45 


X 


X 


V 


V 


V 


7.52 


38.09 


-6.48 


58.17 


9 


EPMl-R 


2.5 X 10" 


0.45 


V 


V 


V 


V 


V 


8.04 


39.59 


-8.23 


58.69 


10 


EPMl 


2.5 X 10-" 


0.45 




X 


V 


V 


V 


7.76 


38.82 


-8.15 


58.41 


11 


PMl 


2.5 X 10"" 


0.45 




X 


X 


V 


v/ 


7.76 


38.82 


-8.15 


58.41 


12 


El / 


2.5 X IQ-" 


0.45 




X 


V 


X 


X 


10.33 


41.29 


-3.62 


60.98 


13 


EPM2-R 


1 X 10-3 


1.8 


V 




V 


V 




7.59 


40.13 


-9.25 


58.85 


14 


EPM3-R 


1 X 10"" 


0.18 


V 


V 


V 


V 


V 


8.24 


40.39 


-7.79 


58.50 


15 


EPM4-R 


5 X lO-'"" 


0.09 




V 


V 


V 


V 


8.29 


40.03 


-7.51 


58.25 


16 


EPM5-R 


2.5 X 10-^ 


0.045 




V 


V 


V 


V 


8.78 


40.33 


-5.77 


58.43 


17 


EPM6-R 


1 X 10--'' 


0.018 


V 


v/ 


V 


V 




9.11 


40.50 


-5.30 


58.37 


18 


EPM7-R 


5 X 10-'' 


0.009 


V 


V 


V 


V 


V 


9.29 


39.97 


-5.14 


58.24 



Note. Model names (except for Model 1, 2, 7 and 8) indicate the activated physics (symbol y/) in the simulations as detailed in Column 
6-9. For example, in EG and El models only mechanical energy feedback is allowed, while in PMO and PMl models momentum and mass 
are considered, but not mechanical energy. Wc adopt 0.2 for the peak EM efficiency cq. 

°77 = 2ewC^/i;^ in v~\g unit where ti„ = 10, OOOkm/s 

''Models with mass, energy and momentum added to the bottom layer. 

'^iglj = iglj opt/^Edd where Lbh is the SMBH luminosity in the optical band after absorption, i.e., as it will be seen from infinity. 

'^AEv, = CwAMbhc^, in erg, where AMbh is in Mq units. 

These models correspond to the models AO, MAO, Al and MAI in Papers I and III, but calculated with some difference in the 
numerical grid spacing. 

^ The model similar to one adopted in SDMH05. 



mimic the com mon treatment of mechanical feedback 
(e.g. SDMH05. iDi Matteo. Soringel fc Hern"^^ [20051 
JNB09). Instead, models indicated with the symbol x in 
Column 5 have a distributed feedback as in Papers I and 
III where we attempt to estimate the gradual deposition 
of mass, energy and momentum taken from the outflow- 
ing wind and going into the ambient g function of 
radius. We then build several models which neglect one 
or two of physical process, i.e., mass, energy, and mo- 
mentum feedback in order to study their effects showing 
the inclusion of each term in Column 6-9. For example, 
Model 3 (EPMO-R, in bold face) distributes the mechan- 
ical feedback only into the bottom layer, and includes 
the radiation feedback and all physical terms, i.e., mass, 
energy, and momentum, in the mechanical feedback. On 
the other hand, Model 6 (EO) adopts a treatment similar 
to that in SDMII05 as it assumes the same mechanical 
feedback efficiency, only includes the mechanical energy 
feedback and distributes it into the bottom layer of the 
grid, neglecting the mass and momentum added back 
into the flow. 

Models 1-6 adopt the standard (high) efficiency ew = 
5 X 10-3, as SDMH05 and JNB09 and Models 7-12 as- 
sume a factor of 20 lower efficiency, perhap s in better 
accord with observation ally based estimates (iMoe et al.l 
120091 lArav et al.l l2010( ) and Models 13-18 adopt other 
efficiencies to show how final properties depend on the 
assumed mechanical efficiency. 



3.1. High efficiency models 

To mimic the common treatment (e.g., SDMII05, 
JNB09), we build the Model 6 (EO) that only includes 
the Mechanical energy feedback with the standard (high) 
efficiency Ew = 5 x 10"'^. In this model, we estimate the 
mass inflowing to the SMBH, convert it to energy with 
the given efficiency, and add this energy only into the 
bottom layers of the surrounding gases. For compar- 
ison. Model EPMO-R has identical efficiency but adds 
also mass and momentum to the bottom layers using 
Equations (5a-d) with 77 — 9, as appropriate for the cho- 
sen efficiency and wind velocity of Uw,io = 1- These two 
models are shown as blue and green lines in Figure [TJ 
We see that allowing for momentum and mass feedback 
reduces the black hole growth by a factor of 1000. The 
more consistent model has a much lower final X-ray lu- 
minosity and final SMBH Eddington ratio. The effect of 
including or not including radiative heating is relatively 
minor, as can be seen by comparing Models 3 and 4 or 
1 and 2. Also the mechanical energy feedback is consid- 
erably less important (as expected) than the momentum 
input, as can be seen by comparing Models 4, 5, and 6. 
Finally it might be thought that some of the effects ob- 
served in these comparisons are due to the change from 
Paper I of adding the feedback to the bottom layers alone 
in the present simulations, rather than over a distributed 
range of radii to mimic the effects of due to a broad line 
wind. But comparison between Models 1 and 3, where 
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Fig. 1. — Models 1-6 with constant and high mechanical efficiency 
e„ = 5 X 10~^ (ri = 9). From top to bottom, the SMBH luminosity, 
X-ray luminosity, mass accreted on the central SMBH, and star 
formation rate are shown with different line types and colors as 
indicated in the third panel. Note how the model that excludes 
momentum feedback, "EO", has by far the highest growth of the 
central SMBH. 

there are only small differences (and Model 1 is identical 
to AO of Paper III) , shows that the differences which may 
be attributed to distributed feedback are small. In sum- 
mary, examination of Models 1-6 shows that including 
momentum drastically increases the effects of feedback. 

3.2. Low efficiency models 

Next, we turn our attention to Models 7-12 which have 
a much lower mechanical efficiency than typically as- 
sumed and it is at a level better in accord with existing 
(and highly imperfect) observational indications. The 
value for the dimensionless parameter rj in these cases is 
only 0.45 (i.e. of order unity), so that we expect that 
inclusion or exclusion of the mass and momentum input 
will make relatively less difference. What do we find? In 
fact, the differences are reduced by about half an order of 
magnitude (0.5 dex), but it remains true that the Model 
12 (like Model 6), without either momentum feedback or 
radiation, has an unacceptably large growth of the cen- 
tral SMBH and an unacceptably large final SMBH lumi- 
nosity, as shown in Figure [5] Models 6 and 12 also show 
thermal X-ray emission greater than 10^^ erg/s, which is 
on the upper side of what is typically observed in normal 



o 

2 






12 



t (Gyr) 



Fig. 2. — Models with constant and low (observation based) me- 
chanical efficiency, ew = 2.5 X 10"'* {rj = 0.45). From top to 
bottom, the SMBH luminosity. X-ray luminosity, mass accreted on 
the central SMBH, and star formation rate are shown, colors and 
line types as in Figure [T] Again momentum feedback is the most 
important physical process in protecting the central SMBH from 
excessive growth. 

elliptical galaxies. 

We summarize the properties of Models 1-12 in Fig- 
ure [3] showing the present-day (14 Gyr) SMBH mass in 
solar mass versus X-ray gas and optical stellar luminosi- 
ties. We show the high efficiency models (Models 1-6) 
in blue and the low efficiency models (Models 6-12) in 
red. The fiducial models, Models EPMO-R and EPMl-R 
with mass and momentum feedback, and Models EO and 
El that only include energy feedback, are marked with 
their model numbers. As discussed above, including the 
momentum and mass feedback not only significantly re- 
duces the SMBH growth but also results in a much lower 
final X-ray luminosity and final SMBH Eddington ratio. 

3.3. Wind efficiency dependence of bottom-layer models 

We test eight different values of Ew for models ranging 
from 5 X 10~^ to 5 x 10~^ for models with the bottom- 
layer treatment and all feedback physics activated (i.e. 
Models EPM#-R). These models correspond to Models 
3, 9, and 13-18 in Table 1. We summarize the results at 
the epoch of 14 Gyr in Figure SI where the least square 
linear fits of several global quantities of interest are also 
given. Note here that the growth of SMBH mass and 
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Fig. 3. — Distribution of all models in the mass-luminosity dia- 
gram measured at z=0, where different colors show different wind 
efficiencies and r] values (tw = 5 X fO~^ and rj = 9 for blue, 
e-w = 2.5 X fO~* and rj = 0.45 for red). The distribution of models 
in the Eddington luminosity - effective SMBH luminosity plane is 
shown in the bottom panel. Five diagonal lines (from top to bot- 
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the BH luminosity Eddington ratio are decreasing func- 
tions of ew If the feedback efficiency is low, too much 
mass is accreted to the central SMBH, as expected. In 
the case of the gas mass and the predicted X-ray lu- 
minosity of the hot ISM, they are decreasing functions 
with increasing wind efficiency but with large scatters. 
We also show the total mechanical feedback energy (i.e. 
AE„ — EwAA/bhc'^) in the last panel, which increases as 
the wind efficiency increases. In this case as in several 
of the others, while the behavior is approximately mono- 
tonic in the direction expected, the approximate power 
law index is much less than unity, since larger efficiency 
gives a larger value for our dimensionless parameter, 77, 
and thus a smaller fraction of the inflowing gas is actually 
accreted onto the central SMBH. 

3.4. The effect of mass removal from the circumnuclear 
disk on purely radiative models 
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Fig. 4. — Dependencies of present-day, global quantities of 
EPM#-R models in Table 1, as a function of mechanical efficiency 
ew From top to bottom, the SMBH mass growth, BH Eddington 
ratio, galaxy gas mass inside 10 Rc, X-raygas luminosity, and total 
wind feedback energy are shown. The linear fits to the data are 
shown in dotted lines, and the fitting results are shown in each 
panel. As expected, the assumption of a higher wind energy effi- 
ciency does correspond to greater feedback effects, but at a much 
less than linear rate. 

In line with the present exploratory discussion, it is 
of some interest also to check the effects of different 
amounts of mass removal from the circumnuclear disk 
via disk wind, in the case of purely radiative models. In 
fact, we recall that in the purely radiative models pre- 
sented in Paper I (such as model RBO in Table 2 therein) 
we do not add mechanical feedback to the equations of 
hydrodynamics, but the mass, momentum and energy 
fluxes of the nuclear wind (and of the jet) are nonetheless 
computed, in order to satisfy Equations [l]-[5] for assigned 
mechanical efficiency and fiducial nuclear wind velocity. 

Therefore, purely radiative models depend indi- 
rectly on the assumed mechanical efficiency, with high- 
efficiency models ejecting a larger fraction of the gas 
from the circumnuclear wind, and therefore reducing the 
amount of gas available for accretion on the SMBH. Here 
we compare the evolution of the purely radiative model 
RBO in Paper I (a model with radiative efficiency O.I 
and with high constant mechanical efficiency 5 x 10^"^), 
with an identical purely radiative model, in which the 
mechanical efficiency has been reduced to zero, therefore 
excluding mass loss from the circumnuclear disk. 

The situation is illustrated in Figure [SJ where the left 
panels refer to model RBO, and the right panels to the 
model without mass ejection from the nuclear disk. In 
the top panels we show the time evolution of the total 
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Fig. 5. — Time evolution of relevant mass budgets (top pan- 
els) and corresponding mass rates (bottom panels) in two purely 
radiative models (without mechanical feedback), differeing in the 
treatment of the circumnuclear disk mass budget. The model on 
the left panels is Model RBO (see Table 2 in Paper I), while the 
model in the right panels is identical in all its properties to RBO, 
except that no mass is lost by the circumnuclear disk. Top panels: 
total mass accreted by the central SMBH (black), of the total mass 
of ISM ejected at 10 Rc (green), and of the total mass in new stars 
accumulated within 10 Rc (red). Bottom panels: the correspond- 
ing mass rates are identified by same colors as in top panels. The 
gas production of the passively evolving stellar population steadely 
declines from ^ IOMq/jt at the beginning down to less than IMq 
at the end. 

mass accreted by the central SMBH (black line), of the 
total ISM mass ejected by the galaxy as a galactic wind 
(green line) , and finally of the accumulated mass in new 
stars (red line). In the bottom line, the corresponding 
rates are shown and identified with the same colors. 

Unsurprisingly, the SMBH grows significantly more (by 
a factor of ~ 2) in the model RBO without nuclear wind 
mass loss (logAMBn/Af© ~ 9.78) than in the model 
with mass ejection (logAMBH/-^© — 9.45). The major 
difference in the accretion history of the two models is 
particularly evident in the first Gyr of evolution, when 
large amounts of gas flow on the central region of the 
galaxy. Note how the SMBH mass of the model with- 
out nuclear mass ejection (right panels) reaches a value 
similar to the SMBH mass of model RBO (left panels) 
at the end of the simulation. As a consequence, the gas 
near the SMBH is gravitationally more bound in the first 
model - especially at early times when the mass losses are 
significant. As can be seen, the star formation history in 
the two models is almost parallel to their SMBH accre- 
tion, and the larger radiative energy output in the model 
without nuclear mass ejection is accompanied by a larger 
starburst at early times, with a final mass of new stars 
of log AM*/M0 ~ 10.5 (red lined), to be compared with 
logAAf,/M0 ~ 10.36 in RBO model without disk mass 
ejection. Finally, consistently with the larger energy in- 
put of the model shown in the right panels, the galactic 
wind expelled a total ISM mass of log AM^/M© ~ 10.4 
in the model without the disk wind, to be compared to 
log AMw/Mq ~ 10.3 in RBO model. 




t (Gyr) 



Fig. 6. — Luminosity versus time for an axisymmetric A model 
with €m = 2.5 X 10~*. Above, the AGN luminosity for half a 
Gyr. Below, the AGN luminosity plotted for a shorter time show- 
ing the highly variable nature of the accretion events in two di- 
mensions. The BH accretion is much more stochastic than the 
one-dimensional case, but the distribution of Eddington ratios is 
quite similar. 

Again, this very simple experiment shows how differ- 
ent treatments in the mass balance equations used to 
describe SMBH can leads to significantly different evolu- 
tionary histories (cf. also iSoker fc Pizzolato 2005). 

4. TWO-DIMENSIONAL MODEL COMPARISON 

One dimensional models continue to be very useful in 
establishing the basic physical processes that are relevant 
for AGN feedback in giant elliptical galaxies. However, 
one dimensional models are not able to capture impor- 
tant properties of the actual systems, including the con- 
vective, Rayleigh- Taylor, and Kelvin-Helmholtz instabil- 
ities. One dimensional models must also rely on a param- 
eterization of the global deposition of mass, energy, and 
momentum via the disk wind, while higher dimensional 
models are able to simulate the evolution of the wind self- 
consistently. We discuss below two-dimensional models 
where we have taken exactly the same galaxy model and 
feedback characteristics to allow comparisons that are 
easy to understand. 

There have been many numerical simulations 
of BH accretion and the subsequent effects on 
the galaxies containing resulting AGN. How- 
ever, efforts to date divide into th r ee categories. 
i Di Matteo. Spri ngel fc Hernguist' (2005"), 'De Buhr et all 
(120091) ■ and [Johansson. Naab fc Burkert . (,2009|) are 
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Fig. 7. — A snapshot from an axisymmetric simulation showing a cold blob falling to the center of the galaxy. On the left, log gas density 
in number of protons per cubic centimeter. In the center, log sound speed in kilometers per second. On the right, the radial velocity in 
kilometers per second. The x and y axes are logarithmic in the distance to the SMBH. The cold gas was produced by enhanced cooling 
in an overdense quasi-spherical shell with a covering fraction of about one third of the sphere. The gas quickly collapses to a ring with a 
small covering fraction and/or fragments as it freely falls to the center of the simulation. 
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Fig. 8. — Black hole mass as a function of time for A models with 
Ew = 10~^. The blue line is the fiducial case with all physics in- 
cluded. The green line has mechanical momentum injection turned 
off so that mechanical feedback is purely via energy; the SMBH 
grows prodigiously, indicating that mechanical momentum feed- 
back is by far the dominant process in limiting the growth of the 
black hole. The red line has mechanical energy injection turned 
off, leaving mass and momentum injection unchanged; the SMBH 
grows somewhat more than the fiducial case, indicating that energy 
feedback plays some role in limiting SMBH growth, albeit a sub- 
dominant one. The cyan line has eem set to zero so that there is 
no energy or momentum feedback due to radiation from the central 
SMBH; this is indistinguishable from the fiducial case, indicating 
that radiative feedback plays essentially no role in limiting SMBH 
growth. 



examples where the simulations cover length scales from 
~ 100 pc to tens of kpc and timescales from a fraction of 
a Myr to several Gyr. Galactic length and timescales are 
resolved, but the BH accretion and feedback processes 



are con sidered to be sub-resolution. iKurosawa fc Progal 
(j2009airbf ) are examples of multi-dimensional simulations 
that cover the length scales from a few AU to ~ 1 pc. 
Length and timescales relevant to BH accretion are re- 
solved, but these simulations do not approach approach 
galactic length o r timescales, and infa ll ra tes are taken 
as giv en. Finallv. iHopkins et al.l (|2009( ) and lLevine et al.l 
(|2008D are examples of a multi-resolution studies of 
BH accretion involving progressively higher spatial 
resolution simulations run for progressively shorted 
times. The highest spatial resolution simulations go 
down to a fraction of a pc, run for about one mega 
year of simulation time. These simulations spatially 
resolve the accretion process, but do not reach galactic 
timescales. Therefore they cannot self-consistently 
calculate the effect of AGN feedback on the gas in the 
galaxy as a whole and subsequent BH accretion. 

The present work is the only attempt of which we are 
aware to simultaneously resolve the inner length scales 
relevant to BH accretion (a few pc), outer length scales 
relevant to galaxies (tens of kpc), inner timescales rele- 
vant to BH accretion (a few years) , and outer timescales 
relevant to galaxies and stellar evolution (10 Gyr). How- 
ever, the region inside of 1 pc including the disk and BH 
itself are still treated as sub-resolution physics. 

A full description of the two-dimensional simulations 
and an analysis of the similarities and differences be- 
tween the one- and two-dimensional models is forth- 
coming. Briefly, we use the Zeus hydrodynamics code 
(jStone fc NormanI 119921 ) in spherical coordinates with 
log-spaced radial bins with Ar/r = 0.1. We have ex- 
tended the code to include appropriate mass, energy, and 
momentum source terms corresponding to stellar evolu- 
tion, star formation, type la and type H supernova feed- 
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back, radiative and mechanical fe edback from AGN ac - 
tivity. See CO07, Papers I and Ill. lSazonov et all (|2005[ ) 
for a full description of the input physics, which are car- 
ried over in all respects except that we have omitted the 
radiation pressure on dust. We require the cells to have 
an aspect ratio of one, giving 30 angular cells. Resolution 
studies have shown little difference in the SMBH accre- 
tion as a function of time as long as the opening angle of 
the disk wind is resolved. 

The major differences between the one-dimensional 
code and the two-dimensional code are in the way that 
the two codes handle angular momentum and the disk 
wind from the AGN. 

The one-dimensional simulations did not permit the 
simulated gas to have nonzero angular momentum. The 
2D simulations assume axisymmetry, but compute the 
velocity in the (j) direction. We must assume an angular 
momentum profile. In the present simulations, we avoid 
forming a rotationally supported gas disk by choosing the 
radius of centrifugal support to be inside the innermost 
grid cell. This allows us to avoid specifying an ad-hoc 
prescription for angular momentum transport. 

The net specific angular momentum of the stars pro- 
viding gas in the simulation is assumed to be: 



1 



1 



R 
j 



d 



(9) 



where R is the distance to the z axis. This parameteriza- 
tion gives solid body rotation at small radii and constant 
specific angular momentum at large radii. The first term 
prevents the rotational velocity from exceeding fa — at 
intermediate radii, there may be a region with constant 
velocity. When gas is created in the simulation by stel- 
lar evolution, it is given this angular momentum profile. 
The subsequent evolution of the gas velocity on the com- 
putational grid is governed the standard fluid dynamics 
conservation laws. 

The one-dimensional code employs a phenomenologi- 
cal model to determine the radius at which energy, mass 
and momentum from the AGN-driven disk wind are de- 
posited in the simulation grid. This model depends on 
an assumed instantaneous jet opening angle. The two- 
dimensional code also requires an assumption about the 
angular dependence of the energy, mass, and momentum 
injected by the disk wind at the edge of the simulation 
grid. Once conserved quantities have entered the sim- 
ulation grid, the two-dimensional code self-consistently 
calculates the time evolution of the material from the 
disk wind; a separate phenomenological model is not re- 
quired. 

For the A models, the opening angle of the jet is chosen 
so that the disk wind covers tt steradians, giving a linear 
opening half-angle of cos~^ | ~ 41°. The opening angle 
does not depend on the BH luminosity in the A models. 
The ID models simply require the jet opening angle as 
a parameter, but the 2D models require that the flux of 
material be fully specifled as a function of angle from the 
z axis. We use 

where q is mass, energy, or radial momentum, £7 is solid 
angle, and 9 is the angle from the z axis. This parameter- 
ization gives half of input material within a half opening 



angle of ~ 41°. 

For the present purpose, the primary result from the 
two-dimensional models is that the qualitative conclu- 
sions already drawn from one-dimensional models re- 
main valid. The dominant physical mechanism regu- 
lating black hole growth is momentum injected by the 
broad-line wind. The energy provided by the mechani- 
cal wind has a noticeable but comparatively small effect. 
The effect of other feedback mechanisms is much smaller 
than either the mechanical momentum or mechanical en- 
ergy- 
Figure [6] shows the AGN luminosity versus time for 
one of the two-dimensional models with a mechanical 
efficiency of ew = 0.001 (corresponding to the one- 
dimensional model EPM2-R). The primary difference in 
the SMBH growth between the one- and two-dimensional 
models is that the two-dimensional models have much 
more stochastic growth. There quiescent periods are not 
as quiescent, and the spacing between the major bursts 
is not as regular in time. Both of these are due to insta- 
bilities present in multiple dimensions: quasi-spherical 
shells of cold gas are able to fragment and fall into the 
center bit by bit rather than as a single large shell. In the 
one-dimensional simulations, bursts of accretion form a 
hot central bubble that is able to prevent further accre- 
tion until the hot bubble cools — this often leads to very 
regular spacing of accretion events in time. 

In two-dimensional simulations, a similar hot bubble 
is formed, but cold gas is able to reach the center via 
Ray leigh- Taylor and convective instabilities. An exam- 
ple of this is shown in Figure [71 Hot gas simply moves 
out of the way leading to much more stochastic SMBH 
accretion with bursts much more closely spaced in time. 

Figure [8] shows the SMBH mass versus time for several 
two-dimensional simulations where each physical process 
is turned off in turn, allowing us to identify which ones 
are negligible and which ones play a dominant role in reg- 
ulating SMBH growth. Without mechanical momentum 
injection, the SMBH grows in a fashion only limited by 
^Edd- Without mechanical energy injection, the SMBH 
grows about a factor of two faster than the fiducial case. 
Mechanical energy plays a role, but it is much less im- 
portant than mechanical momentum input. 
Turning off all radiative feedback processes by setting 
= has little effect on the SMBH growth. Mak- 
ing this choice elimin ates gas h eating as computed by 
the expressions in Sazon ov et al.l (.2005) , momentum pro- 
vided by the absorption of those same photons, as well 
as momentum provided by electron scattering that de- 
termines the Eddington limit. The code does not im- 
pose the Eddington limit — it allows the accretion to be 
limited self-consistently by adding the Eddington force 
to the momentum equation. Therefore setting eem = 
means that the SMBH would not be limited by radiative 
momentum. In spite of this, the mechanical feedback 
is able to keep the accretion rate to physically plausible 
values. The actual optical depth in our simulation for 
electron scattering is typically small compared to unity. 
This is consistent with observations which show that only 
a minority of AGNs are "Compton thick" . 

5. DISCUSSION 

The primary purpose of this paper is to quantitatively 
show, based on one- and two-dimensional computations. 
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exactly which processes are most important during AGN 
feedback episodes; wliicli processes are most useful in 
protecting the central SMBH from excessive mass growth 
and which have most effect on the ambient galaxy. Af- 
ter a central outburst, the mechanical energy must be 
communicated to the ambient gaseous fluid by a wind 
and we, in fact, see these winds in luminous galaxies la- 
beling them "broad-line regions" with outflow velocities 
observed to be ^ 10,000 km/s covering conical regions 
subtending 20-25 % of the sky. These winds must carry 
mass and radial momentum to the ambient fluid, reduc- 
ing thereby the mass deposited on the central SMBH and 
adding a driving component which cannot be reradiated 
away by thermal processes. Equations 5a-d summarize 
the physics, with the dimensionless parameter rj, indicat- 
ing the importance of mass and momentum outflows. 

In the case of an assumed high mechanical efficiency 
(ew = 0.005), we find that, if we suppress the mass 
and momentiim input, then the SMBH grows by over 
a factor of 100 more than if momentum and mass flux 
were properly included in the calculation, and reaches 
masses > lO^^M© in both one- and two-dimensional cal- 
culations. Turning on or off the energy input has rel- 
atively much less effect, altering the SMBH growth by 
roughly a factor of 2. Ignoring the mass and momentum 
feedback inputs also leaves the galaxy with a central op- 
tical luminosity from the AGN which is orders of magni- 
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